      SUBROUTINE SETUP
      RETURN
      END
C   HS71
      BLOCK DATA
      INCLUDE 'BLOCO.INC'
      INTEGER I,J
      DATA NAME/'HS71'/
      DATA (X(I),I=1,4)/1.D0,2*5.D0,1.D0/
      DATA N/ 4/ , NH/1/ , NG/ 9/
      DATA DEL0/1.00D0/ ,TAU0/1.D0/ ,TAU/.1D0/
      DATA (GUNIT(1,J),J=3,10)/8*1/,(GUNIT(3,I),I=3,10)/4*1,4*-1/,
     F     (GUNIT(2,I),I=3,10)/1,2,3,4,
     F     1,2,3,4/
     F     ,((GUNIT(I,J),I=1,3),J=0,2)/-1,0,0,-1,0,0,-1,0,0/
      END
      SUBROUTINE EF(X,FX)
      INCLUDE 'FUCO.INC'
      DOUBLE PRECISION X(NX),FX
      ICF=ICF+1
      FX=X(1)*X(4)*(X(1)+X(2)+X(3))+X(3)
      RETURN
      END
      SUBROUTINE EGRADF(X,GRADF)
      INCLUDE 'FUCO.INC'
      DOUBLE PRECISION X(NX),GRADF(NX)
      ICGF=ICGF+1
      GRADF(1)=X(4)*(X(1)+X(1)+X(2)+X(3))
      GRADF(2)= X(1)*X(4)
      GRADF(3)=1.D0+X(1)*X(4)
      GRADF(4)=X(1)*(X(1)+X(2)+X(3))
      RETURN
      END
      SUBROUTINE EH(I,X,HXI)
      INCLUDE 'FUCO.INC'
      INTEGER I
      DOUBLE PRECISION X(NX),HXI
      CRES(I)=CRES(I)+1
      HXI=X(1)**2+X(2)**2+X(3)**2+X(4)**2-40.D0
      RETURN
      END
      SUBROUTINE EGRADH(I,X,GRADHI)
      INCLUDE 'FUCO.INC'
      INTEGER I,J
      DOUBLE PRECISION X(NX),GRADHI(NX)
      CGRES(I)=CGRES(I)+1
      DO      100       J=1,4
      GRADHI(J)=2.D0*X(J)
  100 CONTINUE
      RETURN
      END
      SUBROUTINE EG(I,X,GXI)
      INCLUDE 'FUCO.INC'
      INTEGER I
      DOUBLE PRECISION X(NX),GXI
      IF ( GUNIT(1,I+NH) .EQ. -1 ) CRES(I+NH)=CRES(I+NH)+1
      IF(I .GT. 1)      GOTO 100
      GXI=X(1)*X(2)*X(3)*X(4) -25.D0
      RETURN
 100  CONTINUE
      IF(I .GT. 5)      GOTO 200
      GXI=X(I-1)-1.D0
      RETURN
  200 CONTINUE
      GXI=-X(I-5)+5.D0
      END
      SUBROUTINE EGRADG(I,X,GRADGI)
      INCLUDE 'FUCO.INC'
      INTEGER I,J
      DOUBLE PRECISION X(NX) ,GRADGI(NX),T1,T2
      IF(I .GT. 1)    GOTO 100
      T1=X(1)*X(2)
      T2=X(3)*X(4)
      GRADGI(1)=T2*X(2)
      GRADGI(2)=T2*X(1)
      GRADGI(3)=T1*X(4)
      GRADGI(4)=T1*X(3)
      RETURN
  100 CONTINUE
      DO      300      J=1,4
      GRADGI(J)=0.D0
  300 CONTINUE
      IF(I .GT. 5)          GRADGI(I-5)=-1.D0
      IF(I .LE. 5)          GRADGI(I-1)=1.D0
      RETURN
      END


